Electrostatic effects in DNA stretching 
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The response of a semiflexible polyelectrolyte chain to stretching in the regimes of moderate and 
weak screening is studied theoretically, with a special focus on DNA experiments. By using the 
nonlinear Poisson-Boltzmann description of electrostatic self-interactions of the chain, we explicitly 
demonstrate the applicability of the concept of effective charge to certain aspects of the problem. 
This charge can be extracted from the far-field asymptotic behavior of the electrostatic potential 
of the fully aligned chain. Surprisingly, in terms of the effective charge, the electrostatically renor- 
malized persistence length can be formally described by the classical Odijk-Skolnick-Fixman (OSF) 
formula, whose domain of applicability is normally limited to linearized Debye Huckel (DH) ap- 
proximation. However, the short scale behavior of the chain in the nonlinear regime deviates from 
the of DH-based result, even upon the charge renormalization. This difference is revealed in the 
calculated stretching curves for strongly charged DNA. These results are in good agreement with 
the recent experiments. In the limit of weak screening we predict the elastic response to have a 
distinctive two-stage character, with a peculiar intermediate " unstretchable" regime. 

PACS: 82.35.Rs, 87.14.Gg 



I. INTRODUCTION 



Single molecular micromechanical experiments are 
among the major tools of modern biophysics. Pioneered 
by the studies of double-stranded DNA (dsDNA)[l], 
these techniques are now applied to a wide range of bi- 
ologically relevant problems [2] -[4]. Theoretically, it has 
been demonstrated that the observed response of dsDNA 
to stretching (i.e. elongation vs. pulling force) can be 
successfully described by semiflexible Worm-Like Chain 
(WLC) model [5] [6], over a considerable range of applied 
forces. Within this model, the molecule is viewed as a 
rigid rod subjected to thermal fluctuations, and the only 
free parameter is its persistence length l p which is pro- 
portional to the bending modulus of the rod. 

Since DNA is strongly charged, the effects of the elec- 
trostatic interactions should in principle play a significant 
role in its elastic response. However, they are typically 
suppressed due to the strong screening under physiologi- 
cal conditions. In addition, it is implicitly assumed that 
the electrostatic effects can be incorporated by renor- 
malization of the parameters of the effective models, e.g. 
persistence length l p in WLC. This is indeed a valid as- 
sumption in certain regimes, as was shown in works of 
Odijk, and Skolnick and Fixman (OSF) [7]. They have 
demonstrated that within linearized Debye-Huckel (DH) 
approximation, there is an additive electrostatic correc- 
tion to the persistence length of a semi-flexible WLC, 
and this term scales quadratically with Debye screening 
length r s . However, the classical result has to be sig- 
nificantly revised in the regime of highly flexible chains, 
which was in focus of numerous studies over the last two 
decades [8]- [11]. 

Another limitation of OSF picture is that it is only 
valid on the length scales larger than r s . When the chain 
is strongly stretched, its properties on shorter scales are 
being probed. Within DH theory, this regime can be 



described in terms scale- dependent electrostatic rigidity 
[12] [6]. Finally, OSF result is known to be inadequate 
outside of DH approximation, i.e. in strongly nonlinear 
regime. Compared to other aspects of the problem, its 
non-linear generalization has attracted only a limited at- 
tention since the early works back in '80s [13] [14]. In this 
paper, we study the electrostatic effects on the elastic re- 
sponse of strongly charged WLC, in the regimes of mod- 
erate and weak screening. In these regimes, DH approx- 
imation may no longer be applicable, and the analysis of 
the non-linear Poisson-Boltzmann equation is necessary. 
While stated as a generic problem, our study is especially 
important in the context of DNA stretching experiments. 

The plan of the paper is as follows. In Section II, 
we formulate the problem and review the results ob- 
tained within DH theory. We revisit the problem of scale- 
dependent rigidity, and demonstrate that in the limit of 
weak screening, the chain exhibits a two-stage response. 
These two stages are characterized by the renormalized 
and bare persistence lengths, respectively, and they are 
separated by a peculiar "unstretchable" regime. In Sec- 
tion III, we study the problem beyond the linearized DH 
approximation. Our central result is that the effective 
rigidity in the non-linear regime can be described by clas- 
sical OSF formula, in terms of effective linear charge 
density determined by the far-field asymptotic of chain's 
electrostatic potential. This is the first explicit demon- 
stration of the applicability of the concept of renormal- 
ized charge in the context of the non-linear problem. 
While the notion of the effective charge is not new, in 
the past it was introduced in a purely heuristic manner. 
As follows from our results, the effective charge approach 
has significant limitations and it is not valid on the length 
scale comparable to r s . As a result, the extension curve 
for strongly charged chains is different from the result 
of linearized DH theory, even upon the charge renormal- 
ization. Finally, in Section IV we discuss applications of 
the theory both to ds- and single stranded (ss-)DNA, and 
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compare our results to the existing experiments. 



II. RESULTS OF LINEARIZED THEORY 

In the presence of monovalent ions of concentration n, 
the electrostatic potential in water is well described by 
Poisson-Boltzmann (PB) equation: 



A* (x) - r s 2 sinh* (x) = -4nl B 



P(x) 



(1) 



Here, reduced potential ^ is expressed in units of kTD/e; 



D is the dielectric constant of water; Ib 



2 /DkT 



7 A is Bjerrum length , and r s = l/^JA-Klsn is Debye 
screening length. For low value of <!/, the above non- 
linear equation can be linearized, which corresponds to 
the well known DH approximation. In particular, one 
can obtain Yukawa-type pair potential acting between 
two point charges e: 



V(r) = l B 
kT r 



exp (—r/r a ) . 



(2) 



The effective Hamiltonian of an electrostatically neu- 
tral WLC stretched by an external mechanical force is 
given by [6]: 



Hn = 



d 2 x 

ds 2 



ds 



ds 



(3) 



Here f is the stretching force applied at s — L. Function 
x (s) determines the chain conformation in space, and it 
is subjected to constrain \ dx/ds\ = 1 (i.e. the chain is not 
extendable) . Within DH approximation, the electrostatic 
self-interactions of uniformly charged chain give rise to 
the following new term in Hamiltonian 



H e l = 



2Dl% 



L L 

II 





exp (-|xQ) -xQ')|/r s ) 
|x( S )-x( S ')| 



dsds' (4) 



Here ae/ls is the linear charge density of the chain. Man- 
ning parameter a has the physical meaning of the number 
of electron charges per Bjerrum length. 

According to OSF theory [7] developed within DH 
approximation, the electrostatic interactions lead to a 
higher effective bending modulus of a semirigid chain. 
This results in renormalization of its persistence length: 



l* p = l p + a 2 r 2 s /4l B 



(5) 



For the case of flexible chain, OSF approach has been 
revised by Khokhlov and Khachaturian, who have in- 
troduced the concept of electrostatic blobs [8]. On the 

1/3 

length scales smaller than blob size, £ e = (l 2 l B /a 2 ) , 
the electrostatic effects may be neglected. On the larger 
scales, a flexible polyelectrolyte can be viewed as a 



stretched chain of blobs, and OSF picture can be recov- 
ered upon proper renormalization of its parameters. For 
a semiflexible chains, we can introduce a characteristic 
length which would play the same role as the size of elec- 
trostatic blobs £ e in flexible case: l e = DkT/(er) 2 = 
Is/a 2 . This length defines the minimal scale at which 
electrostatic effects become relevant, and the semiflexi- 
ble regime itself is determined by condition l e <l p . 

When the chain is strongly stretched, one may keep 
only terms quadratic in the transverse components of 
displacement, u,; = x (s) — f (f • x (s)) /f 2 . After going to 
Fourier representation, u s = L -1 / 2 J2 u q exp (iqs), the 

i 

overall effective Hamiltonian H = Ho + H e i can be writ- 
ten as: 



H[u q ] 
kT 



const + 



<? 2 u 9 u_ 



l p q 2 + Ad (q) + 



±_ 

kT 



(6) 

Here, all the electrostatic effects are collected within the 
term A e ; (q), which has the physical meaning of scale- 
dependent electrostatic tension: 



A eJ (q) 



a 

21b~ 



1 + 



1 



q 2 r 2 



log (1 + <z 2 r 2 ) 



(7) 



This electrostatic term was originally introduced by Bar- 
rat and Joanny [12], and later adopted to the case dsDNA 
in the classical paper of Marko and Siggia [6] . By using 
the Equipartition Theorem, they have calculated the dif- 
ference between the end-to-end distance R and the length 
of the fully-stretched chain , L. 



L-R 



(u 2 ) 



+ 00 



dq 



l p q 2 + A el (q) + f/kT 2tt 



(8) 



The calculated stretching curves are shown in Figure 1 . 
Interesting observations can be made by further analysis 
of the above result in the limit of weak enough screening, 
r s » l e = oT 2 Ib- In this case, the overall mechanical 
response has a striking two-stage character. For large 
forces, / S> kT/Al e we recover well-known result for a 
neutral WLC [6]: 



L-R 



I 




1 dq if 






27 ~ 2V 



In addition, one can identify a new moderate force 
regime, kTl e /Ar 2 s < / < kT/Al e . In this case, Eq. (8) 
has two distinct contributions which correspond to the 
integration over wave numbers q much smaller and much 
larger than r~ , respectively: 




kT 



f(l p + a 2 r 2 /Al B ) 




(10) 
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FIG. 1: Stretching curves calculated for semiflexible polyelec- 
trolytes at various conditions. The solid lines are obtained 
within the linearized DH approximation, while the dashed 
lines correspond to the complete nonlinear PB description of 
the electrostatics. The stretching curve of a neutral chain is 
shown with a dotted line. 

The first term here is similar to the interaction-free 
elastic response, Eq. (9), with the renormalized persis- 
tence length I* = l p + a 2 r 2 s /41,B- This is consistent with 
OSF result, Eq. (5). The second term is independent 
of force, and it indicates that stretching virtually stops 
when kTl e /4r 2 < / <C kT/Al e , while (L- R)/L remains 
finite due to fluctuations on the scales smaller than r s . 
We conclude that both the effective bending rigidity and 
the total chain length should be renormalized due to elec- 
trostatic interactions. In other words, the chain becomes 
nearly " unstretchable" when end-to-end distance reaches 
certain value R* < L (as shown in Figure 1). Electro- 
static effects are relevant only if ar s > y/lslp- When 
this condition holds, function g in Eq. (10) has only 
weak (inverse-logarithmic) dependence on the screening 
length: 

L ~ R * ~ I JH (11) 

L 2a log (ar s / y/l B l P ) ]j 1 p 

The chain remains "unstretchable" until the force be- 
comes strong enough to suppress the small-scale fluctua- 
tions ( / ~ kTlp/r 2 ). After that, the crossover to large 



force regime (9), occurs. 

III. NONLINEAR THEORY AND EFFECTIVE 
CHARGE 

Strictly speaking, the linearization of Eq. (1) is only 
possible for weekly charged chains or very strong screen- 
ing. In a more general case, one has to study the complete 
non-linear equation. In order to apply our results for the 
case of strongly-charged polymers like ds- or ss- DNA, 
one has to go beyond the linearized DH approximation. 
Early attempts of such a generalization go back to early 
80s [13] [14]. At that time, it was shown by the numeri- 
cal studies that effective rigidity clearly deviates from the 
linear OSF theory. However, no alternative model had 
been proposed. 

Various aspects of PB equation in cylindrical geometry 
have been studied for several decades. Among other im- 
portant advances, was the discovery of the counter-ions 
condensation by Manning [17]. It was shown that when 
the linear charge density on a thin rod exceeds certain 
critical value, a finite fraction of the counter-ions gets 
localized within a near region of the rod, thus effectively 
reducing its overall charge. This phenomenon occurs at 
critical Manning parameter a c ~ 1. Upon the condensa- 
tion, the residual effective charge of the rod saturates at 
its critical value, <x c /Ib — l/'s- For a strongly-charged 
cylinder with a finite radius vq , a similar phenomenon oc- 
curs with the effective charge becoming a function of as- 
pect ratio, ro/r s . Following Oosawa [18], the idea of sep- 
arating condensed and free counter-ion components, has 
been widely used in literature. Nevertheless, it is hard 
to justify its applicability for our problem. First, there 
is an ambiguity in defining the "near" and "far" region 
of the cylinder with a finite aspect ratio, and therefore 
the effective charge is not well defined. Furthermore, the 
electrostatic interactions are strongest on lengths scales 
shorter than r s , exactly were DH approximation is invalid 
(for a > 1). 

The leading corrections to the free energy due to 
small fluctuation of the chain can be evaluated with- 
out any model assumptions. Namely, one can lin- 
earize PB equation near its solution for a static cylin- 
der (\I>( ) (x)), and couple the perturbations of potential 
4> (x) = 'F (x) — , F(°) (x), to the conformational fluctua- 
tions u s . The displacements of the charged cylinder can 
be described by the local perturbations to the in-plane 
electrostatic moments. In particular, the monopole, 
dipole and quadrupole moments have the following forms, 
respectively: 

ae ( du s \ 2 

Ts = 2T B \-d7) ' (12) 

d s = — u s ; (13) 

IB 



4 



<-B 



o 



ds 



■ (14) 



Here, a = uq + (ro/2r s ) 2 sinh (ro) takes into account 
the fact that the cylinder of radius r with linear density 
density eao/ls, fluctuates in oppositely-charged environ- 
ment with local charge density ne sinh \E , (°- ) (ro) . We have 
limited our consideration to the leading three electro- 
static moments because they are the only ones containing 
displacement u s in orders no greater than 2. Thus, the 
general expression for the correction to the electrostatic 
Hamiltonian (second-order in u) is, 



5H, 



6H e 
5t s 



H. 



(0) 



-T s - 



2L 



ds 



ds+ 



(15) 



Thus, 

5J T^ = E[t7 M - ^ (ro) ) u « u « + Cq2 ^ u - 

(19) 

Here the we have used the fact that the first two terms 
in Eq. (15) are local in du s /ds: 

const ■ / (du s /ds) 2 ds = Cq 2 u q \i- q 

J q 

We can determine coefficient C by noting that the low- 
q behavior of the electrostatic correction to the Hamil- 
tonian, H e i [u s ] should still be dominated by bending 
modes, q 4 u q u_ q , and therefore all the terms contain- 
ing q 2 u q U- q must cancel in that limit: 



L L 



1 

+ 2 



S 2 H et 
<5d s <5d s / 



o o 



d,d,/dsds' 



Here, the first term accounts for the fluctuations of the 
local linear charge density, t s . By symmetry reason, a 
similar term containing quadrupole moment Q s is ab- 
sent (SH e i/5Cl s = 0). The second term in the above 
expression reflects the fact that the cylinder is no longer 
perfectly straight and its end-to-end distance is reduced 

L 

by amount J (du s /ds) ds/2. The last term represents 
o 

the dipole-dipole interactions, and this is the only non- 
local term in the effective Hamiltonian: 



L L 

1 f f ( 8 2 H el 



2 J J \Sd s Sd s 





d s /d S 'dsds' = 



(16) 



& 2 I d 2 ^ (r ) 
1 l B \ dq 2 



(20) 



9=0, 



After substitution of this value of C into Eq. (19), 
we obtain a nonlinear generalization of the electrostatic 
tension term A e ; (q) which enters the earlier results, Eqs. 
(6) and (8): 



A 



(nonlin) 
el 



4<5 2 

(9) = -^- 



4 X) (ro) - ^ (r ) d 2 ^ (r ) 



2dq 2 



9=0/ 



(21) 



By solving Eqs. (1) and (17) numerically, we found 



function A 



(nonlin) 

•X 



(q) for various regimes. Its low-g 



asymptotic may significantly deviate from classical OSF 
result, which is consistent with the earlier studies [13] 
[14]. One could interpret this deviation in terms of renor- 
malized effective charge, by defining the effective Man- 
ning parameter oS e f so that in the limit qr s <C 1, 



L L 



~D J J ^ (r °' S ~ S>) (ds ~ ds ' )2 dsdS ' = 



o o 



= |E(4 1} M-^o 1} (to)) d,d_ 9 



Here we have switch to Fourier representation, and ex- 
pressed the result in terms of the dipolar field <$P (r). 
The latter satisfies the inhomogeneous linearized PB 
equation: 



d 2 + -d r - (q 2 + 1 + 1 cosh^W (r) 



W = o, 



subject to boundary condition: 



r 



1. 



(17) 
(18) 



A 



(nonlin) 
el 



(q) ~ 



{a^f)r s qy 



(22) 



Surprisingly, this effective effective charge appear to have 
a clear physical meaning. With an excellent accuracy, 
a (e//) h as t ne game value as effective Manning parameter 
exoo which can be extracted from the far-field behavior of 
potential *W (r), i. e. 

tf<°> (r) ~ 2 aoo Ko(r/r s ) 

for r > r s [16]. Here K is the modified Bessel func- 
tion which is the solution to DH equation in cylindric 
geometry. 

Figure 2(a) demonstrates a striking correlation be- 
tween oS e $^ and a^, in a wide range of parame- 
ters. Each curve there corresponds to the fixed value 
of aoo. The solid and dashed lines represent a^ e ^^/a 
and aoo/ao respectively, as functions of the aspect ra- 
tio ro/r s . The close agreement between the two types 
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FIG. 2: (a) Correlation between the effective Manning pa- 
rameter a?' e -f-f' obtained from the renormalized rigidity of the 
strongly charged chain, and the one extracted from the far- 
field behavior of the electrostatic potential, a„. Solid and 
dashed line represent / qq , and ao , respectively, 

for the fixed value of a x and variable aspect ratio ro /r s . (b) 
a (eff) ag a f unc tion of aspect ratio for various values of bare 
charge: ao —* oo ( the saturation curve, solid line), ao = 4 
(dsDNA, dashed line) and ao = 1 (ssDNA, dotted line). Note 
that eS ej ^ clearly devites from the result obtained through 
Oosawa-type calculation, Ref. [19]. (c) Effective vs. bare 
charge dependence (solid lines), compared to similar result of 
Oosawa-Manning model (dashed lines). The two sets of plots 
correspond to the limit of an infinitely thin rod (ro = 0), and 
a finite aspect ratio (ro/r s = 0.1), respectively. 
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FIG. 3: A e i (q) obtained from the complete PB equation, for 
various aspect ratios (solid lines). An introduction of the ef- 
fective Manning parameter results in matching of these curves 
with DH (dashed), and OSF (dotted) results, for low q. For 
qr s > 1, the DH approximation underestimates the effect 
of electrostatic interactions (compare the zero aspect ratio 
curves). Both PB and DH curves significantly deviate from 
OSF parabola, indicating that the effective rigidity descrip- 
tion is not valid on length scales smaller than r s . 



of data shows that a^ e ^^ and nearly coincide for an 
arbitrary value of the bare linear charge, as long as the 
aspect ratio ro/r s remains small. Figure 2(b) illustrates 
the dependence of a^ e ^^ on the aspect ratio for the val- 
ues of the bare charge close to those of dsDNA and ss- 
DNA (ao = 4 and ao = 1, respectively). One can see 
that for the practical purposes, a^f^ of ssDNA is close 
to .7, while its value for dsDNA follows the saturation 
curve (corresponding to infinitely-charged cylinder with 
the same radius). Note that the far-field definition of 
the effective charge is unambiguous, and it does deviate 
from Oosawa-Manning result, even for zero aspect ratio 
(see Figure 2(c)) . Furthermore, one can modify Oosawa- 
Manning approach to evaluate the effective charge as a 
function of ro/r s , as demonstrated in Ref. [19]. The re- 
sult of that calculation is also shown in Figures 2(b)-(c), 
and it clearly deviates from our far-field based definition 
ofa< e// ). 

We have explicitly demonstrated that OSF result is 
formally applicable even outside of DH approximation, 
given that the effective charge is properly defined. While 
similar effective charge descriptions have been previously 
utilized in several works [16] [15], they were all introduced 
as purely empirical assumption. In fact, our result is 
rather non-trivial since the electrostatic interactions are 
strongest in the vicinity of the chain, where the effec- 
tive description is not valid. The effective charge picture 
works only because the bending of the chain affects the 
far region much stronger than the near one. Indeed, q 2 
term in Eq. (17) is most relevant for r > r s . 

Furthermore, for q > l/r s the effective charge descrip- 
tion breaks down, and deviations from DH approxima- 
tion become considerable, as shown in Figure 3. Qual- 
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itatively, this can be interpreted as an increase of the 
effective charge with q: as one goes to smaller length 
scales, a( e f^> should shift from its far-field value ctoo, to- 
wards the bare Manning parameter ao- The change of the 
overall shape of A e ; (q) gives rise to modification of the 
stretching curves, compared to the results of DH-based 
theory, as shown in Figure 1. Naturally, these corrections 
are strongest for large ratios r s /l p , which corresponds to 
wider dynamic range of relevant wave vectors. 



IV. DISCUSSION 

The above results can be applied to the experiments 
on both ds- an ssDNA. The two cases are character- 
ized by rather different values of the physical param- 
eters. In particular, the persistence length of dsDNA, 

l p ds ' — 50nm ~ 701 b significantly exceeds that of ssDNA 
chain, l p ss ^ ~ lnm ~ Ib- In addition, their linear charge 
density differ approximately by factor of 4. Finally, while 
WLC model is well established as a standard description 
of dsDNA, the situation with elasticity of ssDNA is not as 
settled: the extensible version of WLC is among several 
proposed models all of which are in reasonable agreement 
with the existing experimental data [20]. Plots in Fig- 
ures 1(a) and (b) roughly correspond to the parameters 
of ds- and ssDNA, respectively. As one can see from 
the Figure, the two-stage behavior and the intermedi- 
ate "unstretchable" regime becomes pronounced when r s 
substantially exceeds the persistence length. Note how- 
ever that this regime is unlikely to be accessible in the 
case of dsDNA, since the double helix would become un- 
stable at such a week screening due to the electrostatic 
repulsion between the strands. For ssDNA, the two-stage 
stretching must be observable at r s > lOOnm. 

While Figure 1 gives only a general idea about the 
mechanical response of semiflexible chain with electro- 
static interactions, a more detailed comparison of our 
results with experiments is presented in Figure 4. The 
plot shows the stretching curves of ssDNA taken at vari- 
ous screening conditions, as reported in fiefs. [2] and [4]. 
These experimental data are fitted by our theory with the 
Debye screening length calculated for each salt concen- 
tration. As a result, the only fitting parameter was the 
bare persistence length l p , and the best fit is achieved at 
l p = 5.5A This value is somewhat smaller that the ear- 
lier estimates of the effective persistence length I* ~ 8A, 
which difference is due to our explicit account for the elec- 
trostatic corrections to chain rigidity. Unlike the plots in 
Figure 1 which for the sake of simplicity were calculated 
by assuming ro/r s = 0, in Figure 4 the non zero diam- 
eter of ssDNA chain is taken into account. Based on 
the known structure of the molecule, we have assumed 
ro = 5A, and this parameter was not used for fitting. 

Even though the theoretical stretching curves were cal- 
culated with WLC model, our approach can be also 
applied to alternative descriptions of ssDNA, such as 




FIG. 4: Comparison of the theoretical results to experimental 
data on elastic response of ssDNA at various screening con- 
ditions. Data points are taken from Refs. [2] (2mM NaCl, 
shown by circles), and [4] (ImM and lOmM PB data are rep- 
resented by squares and diamonds, respectively). The theo- 
retical curves were obtained by taking into account the calcu- 
lated values of Debye length for each case ( 6.9nm, 6nm and 
1.9nm respectively). ssDNA was approximated by a cylinder 
of radius ro = 0.5nm. The only fitting parameter for all three 
curves is bare persistence length l p , whose optimal value is 
close to 5.5nm. 

Extensible Freely- Joint Chain model (EFJC). However, 
this modification is unlikely to produce any significant 
change since within the force range of interest, the com- 
peting models (WLC, EFJC, and discreet version of 
WLC) are known to give nearly identical fitting to the 
electrostatics-free stretching behavior of ssDNA [20] . In- 
directly, this model independence is supported by the 
fact that the above experimental data were also success- 
fully fitted with numerical simulations based on EFJC 
model [4]. In that work, in addition to electrostatic inter- 
actions described within linear DH approximation with 
renormalized charges, the effects of excluded volume and 
hairpin formation were included [15]. Since these addi- 
tional effects are only important for relative small chain 
extension, they were neglected in the context of our work. 



V. CONCLUSIONS 

In conclusion, we have studied the effects of electro- 
static self-interactions on elastic properties of strongly- 
stretched WLC, with application to ds- and ss-DNA at 
moderate and weak screening conditions. By revisiting 
the earlier studied problem of scale dependent electro- 
static rigidity we have demonstrated that in the limit of 
weak screening, the mechanical response of the chain has 
a distinctive two-stage character. The two modes are 
characterized by the bare and the renormalized persis- 
tence length, respectively, and they separated by inter- 
mediate "unstretchable" regime. 
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The central theme of our work is the effect of non- 
linearity of PB equation on the classical results obtained 
within DH approximation. Consistently with previous 
studies of the problem [13] [14], the OSF result for the 
electrostatic renormalization of the persistence length is 
not applicable in the nonlinear regime, and in particu- 
lar the scaling I* — l p <~ r 2 s is not valid. However, we 
have explicitly demonstrated that the OSF formula is 
formally correct if the actual line charge density is re- 
placed with the effective one. The latter is determined 
from the far-field asymptotic of the electrostatic poten- 
tial of the strongly stretched chain. In the light of this 
result, the deviation from OSF scaling is due to the non- 
trivial dependence of the effective charge on the aspect 
ratio ro/r s . 

The renormalized OSF description is however only par- 
tially adequate since the electrostatic corrections to the 



chain rigidity do not follow the result of DH theory on 
scales shorter that r s , even upon the charge renormal- 
ization. This deviation is due to the fact that the effec- 
tive charge itself is a scale-dependent concept. Our solu- 
tion to the complete nonlinear problem is shown to give 
an adequate description to the experimentally observed 
stretching behavior of ssDNA at various salt conditions. 
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